Geometric Phase, Curvature, and Extrapotentials in Constrained Quantum Systems 
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We derive an effective Hamiltonian for a quantum system constrained to a submanifold (the 
constraint manifold) of configuration space (the ambient space) by an infinite restoring force. We 
pay special attention to how this Hamiltonian depends on quantities which are external to the 
constraint manifold, such as the external curvature of the constraint manifold, the (Riemannian) 
curvature of the ambient space, and the constraining potential. In particular, we find the remarkable 
fact that the twisting of the constraining potential appears as a gauge potential in the constrained 
Hamiltonian. This gauge potential is an example of geometric phase, closely related to that originally 
discussed by Berry. The constrained Hamiltonian also contains an effective potential depending on 
the external curvature of the constraint manifold, the curvature of the ambient space, and the 
twisting of the constraining potential. The general nature of our analysis allows applications to a 
wide variety of problems, such as rigid molecules, the evolution of molecular systems along reaction 
paths, and quantum strip waveguides. 
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I. INTRODUCTION 

> 

^\ , We derive an effective Hamiltonian for a quantum system subject to an infinite restoring force. Though our results 
■ are quite general, we are motivated by several specific applications, such as stiff molecular bonds in rigid molecules 
and clusters of rigid molecules, molecular systems evolving along reaction paths, and electrons confined to quantum 
strip waveguides. 

For comparison, consider first a classical system. We have in mind a system which initially occupies any position 
in the full configuration space (called the ambient space) but is subsequently confined to a submanifold (called the 
constraint manifold) by the introduction of a restoring force, which in a certain limit becomes infinite. Here, the 
Hamiltonian is simply the kinetic energy plus the constraining potential, which we assume is constant along the 
constraint manifold. Assuming the initial velocity is tangent to the constraint manifold, it is well known that the 
system's trajectory remains on the constraint manifold and that its motion is determined solely by the form of the 
kinetic energy tangent to the manifold Jl],^| . This kinetic energy, in turn, depends only on the Riemannian metric of 
the constraint manifold. Thus, the motion of the constrained classical system depends only on the internal metric of 
the constraint manifold and is independent of the ambient space, the embedding of the constraint manifold within the 
ambient space, or the details of the constraining potential. It is a remarkable fact, then, that for a quantum system 
this is no longer true. The constrained quantum Hamiltonian depends on the curvature of the ambient space, the 
external curvatures of the constraint manifold, and on properties of the constraining potential. 

It has already been known for some time that a constrained quantum system "senses" the local neighborhood of the 
constraint manifold . As a simple example, consider a quantum system whose motion transverse to the constraint 

manifold is in the ground state. Due to the conservation of the transverse action, the constrained quantum system 
sees the transverse zero point energy as an effective potential on the constraint manifold. (We call this the adiabatic 
potential.) The adiabatic potential also appears classically if the initial velocity of the system has a nonzero component 
normal to the constraint manifold. Classically, the adiabatic potential can always be eliminated by choosing an initial 
velocity tangent to the constraint manifold. Quantum mechanically, however, due to the Heisenberg principle, the 
transverse action and hence the adiabatic potential can never be eliminated. 

The present paper focuses on effects of the ambient space and constraining potential other than the adiabatic 
potential. Following da Costa [p|-|s|,[io| , we assume that the constraining potential has the same form at each point of the 
constraint manifold. The adiabatic potential is therefore constant along the constraint manifold and can subsequently 



be ignored. (In Sect. VIC, we discuss briefly how a small amount of variation in the adiabatic potential can be 
accommodated.) In two noteworthy papers 1^,0, da Costa, using this assumption, derived the effective Hamiltonian 
for a system of n constrained point particles. This Hamiltonian contains two terms. The first is proportional to the 
Laplacian on the constraint manifold, and therefore depends only on the internal metric of the constraint manifold. 
The second, however, is an effective potential, called the extrapotential, which depends not only on the internal 
curvature, but also the external curvature of the constraint manifold. This extrapotential is of order ti 2 and therefore 
vanishes in the classical (and semi-classical) limit. As a simple, yet illustrative, example, consider a system defined on 
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M 3 constrained to lie on a curve. For this system, the extrapotential is — h 2 / (8p 2 ), where p is the radius of curvature. 
This result was obtained by da Costa ||; the same result was obtained earlier by Marcus ||] and Switkes et al. (|] 
for curves in M 2 . Others have also studied this order H 2 extrapotential, including Jensen and Koppe and Kaplan, 
Maitra, and Heller |rjf| . Since the extrapotential depends on the external curvature, it can never be derived from 
a procedure which quantizes the constrained classical system (which depends only on the internal curvature), an 
approach common in the literature of constrained quantum mechanics. (See, for example, the review of DeWitt |L2|.) 

As mentioned above, once the constraining potential is defined at one point of the constraint manifold, the con- 
straining potential at all other points must have the same form. However, this requirement does not completely 
determine the constraining potential at all points since the orientation of the potential is left unspecified. In other 
words, the equipotentials surrounding the constraint manifold can twist in some unspecified manner as one moves 
along the manifold. 

Da Costa fixed the twisting ambiguity by imposing what we call a "no twist" condition on the potential. Physically, 
this condition requires the restoring forces in the neighborhood of the constraint manifold to be normal to the manifold. 
It can be viewed as an extension of the fact that in classical mechanics nondissipative forces are normal to the constraint 
manifold at the point of the manifold itself. Da Costa astutely realized that if the no twist condition were violated, 
the motion transverse to the constraint manifold would be coupled to the motion tangent to the manifold and the 
Schrodinger equation would not separate. 

For some submanifolds there exist no potentials which satisfy the no twist condition. In Ref. Q| da Costa derived a 
local geometric criterion on the external curvature of a submanifold which was necessary and sufficient to determine 
the existence of a non-twisting potential. Unfortunately, several common examples of constrained systems do not 
satisfy this criterion. For example, consider a model of a polymer by n > 2 point particles where the distances 
between each particle i and its neighbor i + 1 are fixed. (This model also applies to the double pendulum.) These 
systems all fail the criterion |13| as does a system of n > 2 point particles constrained to form a rigid body. Even if a 
given submanifold can have a non-twisting constraining potential, there is no guarantee that this potential is the one 
dictated by the physics of the system under consideration. 

The principal objective of this paper is to derive an effective Hamiltonian for a constrained quantum system with 
arbitrary twisting of the potential. The presence of the potential twist leads to several qualitative changes in the 
Hamiltonian. First, the Hamiltonian is no longer a scalar operator, but rather a k x k matrix of operators acting 
on a fc-dimensional vector-valued wave function defined over the constraint manifold. Here, k is the degeneracy of 
the transverse modes with each component of the vector wave function representing a different transverse mode. Of 
course, if one chooses a nondegenerate transverse mode, the Hamiltonian reduces to a single component. 

Perhaps the most interesting consequence of dropping the no twist condition is the emergence of a U (k) gauge 
potential, or connection, in the constrained Hamiltonian. This gauge potential is a coupling between the twisting of 
the potential and the generalized angular momentum of the transverse modes. Modes with no such angular momentum 
are unaffected by the potential twist. This gauge potential is an example of geometric phase and is closely related to 
the phase originally introduced by Berry in the context of adiabatic transport of quantum states (TJ] • It is interesting 
to note that the gauge potential is of order h° and therefore, like the adiabatic potential, is essentially of classical 
origin. 

In addition to creating the gauge potential, the potential twist adds a term to the extrapotential. Unlike the 
extrapotential terms derived by da Costa, the potential twist term is not a scalar function, but a k x k matrix of 
such functions with possible off-diagonal terms coupling the degenerate transverse modes. The potential twist term 
depends on the standard deviation of the angular momentum of the transverse modes and thus disappears for angular 
momentum eigenstates. 

In some applications, the ambient space may not be flat. For example, the internal space of a molecule with n > 2 
atoms is not flat [ jl5| . Constraining such a molecule to a reaction path therefore requires an analysis accounting for 
the ambient curvature. We therefore do not assume in this paper that the ambient space is flat. The effects of the 
ambient curvature are most notable as additional terms in the extrapotential, although it also modifies the curvature 
of the gauge potential. 

This paper has the following organization. In Sect, [nj we introduce many of the key concepts by a simple example: 
that of a system on R 3 constrained to a curve. Section || is purely expository, containing no derivations. Section III 
briefly introduces the general problem. In Sect. IV we focus on the constraining potential. We take care to precisely 
define the requirement that it have the same form at all points of the constraint manifold. We also define a tensor 
which measures the twisting of the potential. In Sect. IV C we specify how the potential is to scale in e, where e — > 
represents an infinite constraining force. The main computations of the paper are in Sect. [v| in which we expand the 



kinetic energy in e and arrive at a preliminary expression, Eq. (5.31), for the constrained kinetic energy. Section VE 



devoted to deriving various expressions for the extrapotential. In Sect. VI, we apply first order perturbation theory to 
transform to a set of degen erate transverse modes, thereby obtaining Eqs. (3.21) - ( |6.23 ) , wh ich are the main results 
of the paper. Section VI C briefly discusses nonconstant constraining potentials. In Sect. VII, we study the geometric 
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origins of the gauge potential and various related connections. We also compute their curvatures. Section VIII contains 



some special cases, including constraint manifolds of codimensions one and two, rotationally invariant constraining 



potentials, and harmonic constraining potentials. In Sect. VIII E, we show that the gauge potential vanishes for 
certain systems with reflection symmetry. Conclusions are in Sect. IX. There are three Appendices. Appendix A 
contains a very brief review of curves in M 3 . Appendix |b] is a review of the second fundamental form. Appendix C 
summarizes an expression we will need for the quantum kinetic energy. 

II. A SIMPLE EXAMPLE: A CURVE IN R 3 

The ultimate objective of this paper is to constrain a quantum wave function, defined on an arbitrary manifold (the 
ambient space), to a (locally) arbitrary submanifold (the constraint manifold) via some general constraining potential. 
Before solving the full problem, however, it is instructive to consider a simple (though certainly non-trivial) , c oncrete 



example of the constraining procedure. We present no derivations here; our results will be justified later in Sect. VIII _B . 

We consider a curve embedded in flat three-dimensional space R 3 and parameterized by its arclength a. (See Fig. |l|.) 
Such a curve is characterized by its curvature n and torsion r. (See Appendix |A].) We take this curve to be the axis 
of a quantum waveguide. That is, there is a tube enclosing the curve such that the potential is zero inside the tube 
and infinite outside. We assume the cross-section of the tube is constant along the curve. More precisely, if we cut 
the tube along a plane normal to the curve (called hereafter a normal plane), the cross-sectional shape of the tube 
is independent of where along the curve we cut. Two such cross-sections have the same shape if one can be rotated 
into the other. This rotational freedom permits the cross-sectional shape to twist as one moves along the curve, even 
if the curve itself is straight. The orientation of the cross-sectional shape is specified by two orthonormal vectors Ei 
and E2 chosen from each normal plane along the curve. The choice of orthonormal frame (Ei,E2) is such that the 
cross-sectional shape (with respect to this frame) is independent of a. In Fig. [|, the cross-section is a triangle with 
no re flection symmetry. Such symmetry is nongeneric and can cause certain terms discussed below to vanish. (See 



Sect. [VIII n ) 

We assume that the transverse dimensions of the tube are small compared to the radius of curvature p = k^ 1 
and the inverse torsion r -1 . We can then separate out the "fast" transverse degrees of freedom and obtain an 
effective one-dimensional Hamiltonian in the "slow" longitudinal, or tangential, coordinate a. To accomplish this 
separation, we pick a transverse mode ^(u 1 ,^, 2 ) of the waveguide. Here (u 1 ,^ 2 ) are the Cartesian coordinates in the 
normal plane taken with respect to the frame (Ei,E2); the quantities (u 1 ,^ 2 ,^) thus coordinatize the tube. The 
transverse mode x( u \ u 2 ) is a normalized eigenfunction of the transverse Hamiltonian H±_ = (ti\ + ir 2 ) / '2 + V±(u , u 2 ), 
where ttj — —ihd/du^, j — 1,2, and V±(u 1 ,u 2 ) is the potential energy which defines the tube. The eigenvalue of 
H±_ corresponding to \ is called the transverse energy. For simplicity, we assume that the transverse energy is 
nondegenerate . 

To lowest order in the width of the tube, an eigenfunction if} of the wave guide has the form 

TP(u\u 2 ,a)= X (u\u 2 ) ( f>(a). (2.1) 

As we take the limit where the transverse dimensions of the waveguide shrink to zero (keeping the quantum numbers 
of the transverse mode fixed), the transverse energy obviously tends toward infinity. However, due to the constancy of 
the cross-sectional shape, this transverse energy, though very large, is itself constant along the curve. We thus subtract 
it off, leaving a residual Hamiltonian Hu , which we call the constrained Hamiltonian. The constrained Hamiltonian 
acts only on 0, resulting in the Schrodinger equation 

H\\4> = E\\4>. (2.2) 

The principal objective of this paper is to determine the form of this constrained Hamiltonian. 

As we will show later, the constrained Hamiltonian is not simply tt 2 /2 where 7T|| = —ihd/da. Rather, there are 
effects from the curvature k and the rate at which the cross-sectional shape twists along the curve. To make this 
latter concept more precise, we define the potential twist S = Ei ■ {dEi2/da) which measures the rotation rate of 
the cross-sectional shape. The potential twist admits the following description. Let denote the angle between the 
principal normal n (see Appendix |X|) and the frame (Ei,E2), specifically, n • Ei = cos 6*, n • E2 = —sin 9. Let 
u) = d9/da denote the rotation rate of the frame (Ei, E2) with respect to n. Then S is related to u> and the torsion 
t by —S = t + u>. Taking S — 0, we obtain the case considered by da Costa in Ref. |J. We next define an angular 
momentum operator A in the tangential direction by A = (u %2 ~ u 2 tti)/2. The constrained Hamiltonian is then 

Hy i =K Vi + V ex , (2.3) 
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where 



K }] =-(7T|,+ 25(A)) 2 , (2.4) 

K,=-|^ 2 + 2S 2 ((A 2 )-(A) 2 ), (2.5) 

and where the bracket notation ( ) denotes the expectation value with respect to the transverse mode x- 

Observe that the tangential kinetic energy K\\ departs from 7r|/2 due to the inclusion of the term 25(A), which 
couples the angular momentum o f th e transverse mode to the rate of potential twist. This term is a gauge coupling, 
a fact we explore further in Sect. I VIll. F or now, we simply note that because the curve is one-dimensional, the gauge 
coupling can be removed from Eq. by a gauge transformation. In the present context, a gauge transformation 

consists of changing the phase of the wave function <f>. This transformation is not without its consequences, however, 
as it will obviously change the boundary conditions which <j> must satisfy. Also, we stress that if the constraint 
manifold has dimension greater than one, it will not in general be possible to remove the gauge coupling by a gauge 
transformation. 

As a final observation on Ku, notice that the gauge coupling is of order H , which indicates that it is essentially a 
classical quantity. This coupling should therefore appear in a classical theory of constraints which takes into account 
the potential twist. 

Turning to the quantity V ex , we note that it is a real- valued function of a, containing no differential operators. For 
this reason, we call V ex an extrapotential. The extrapotential contains two terms, — h 2 K 2 /8 and 2S' 2 ((A 2 ) — (A) 2 ). 
The first of these was derived by da Costa for the case S = |i| . It has the physical effect of attracting (f> to regions 
of high curvature, a fact which may produce curvature-induced bound states in the waveguide. Such bound states are 
of current interest [p^^8[ and are reviewed by Duclos and Exner pjj . The term —h 2 K 2 /S is of order h 2 and therefore 
disappears in the classical (and semi-classical) limit. The second term of V ex , like the gauge coupling in Ku, depends 
on both the potential twist S and the angular momentum A. Notice, however, that it is the standard deviation of 
the angular momentum which appears in V ex . This means, for example, that the second term of V ex vanishes for 
transverse modes which are angular momentum eigenstates. It is interesting to observe that, in contrast to the first 
term, the second term of V ex has the physical effect of expelling the wave function (f> from regions of high twist S. 
Also, the second term is of order h°, which means that, like the gauge coupling in Ku, it survives the classical limit. 



III. INTRODUCTION TO THE GENERAL PROBLEM 

We describe here how the setup in Sect. |0| is modified for the general problem. First, the ambient space in Sect. H 
was assumed to be R 3 . In the general problem, we allow the ambient space to be an arbitrary Riemannian manifold, 
which we denote by A. The kinetic energy of the wave function ip, defined over A, is given by K = — fi 2 A/2, where 
A is the Laplacian on A. Unlike Sect. [Il| the ambient space is not assumed to be flat, and, as we will discover, the 
curvature of the ambient space produces additional terms in V ex . 

Next, we constrain the wave function to lie in the vicinity of a (locally) arbitrary (embedded) submanifold C of A 
with dimension m and codimension d. We call C the constraint manifold. In Sect, [H], the constraint manifold was a 
one-dimensional curve. The curvature and torsion of this curve played a critical role in the analysis. The appropriate 
generalization of the curvature and torsion is the second fundamental form T, which is a rank three tensor. (See 
Appendix 0.) 

In Sect, jy, the constraint was imposed by a hard-wall potential that was infinite outside of a tube and zero inside. 
We then took the limit in which the width of the tube went to 0. In the general problem, we impose the constraint 
by an arbitrary potential V±, subject to a few reasonable conditions. This potential is defined on a set of coordinates 
transverse to the constraint manifold and, for this reason, is called the transverse (or constraining) potential. The 
transverse potential depends on a scaling parameter e which is analogous to the width parameter of the tube; the 
constraint is imposed by taking the limit e goes to 0. 

One of the conditions we do still require of V± is that it be independent of the location on the constraint manifold. 
This condition, as well as a few other minor conditions, are explained fully in the next section. 



IV. THE TRANSVERSE POTENTIAL 
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A. Constancy of the Transverse Potential 



In Sect. we denned the constraining potential by first specifying the form of the potential on a plane normal 
to the curve and then specifying the orientation of this potential at all points along the curve. For the general case, 
we use the same fundamental idea except that now, due to the curvature of the ambient space, we must take care to 
define how we generalize the concept of the normal plane. 

It is useful to consider two separate but related spaces for a given point q on the constraint manifold. The first is 
the linear space of all vectors normal to the constraint manifold. We call this the normal space at q and denote it 
by N q . The second space of interest is the submanifold of the ambient space formed by geodesies emanating from q 
normal to the constraint manifold. We call this the transverse space at q and denote it by U q . The spaces N q and lA q 
are related by the exponential map which takes a vector v £ N q into the point exp v G A. The point exp vei lies on 
the geodesic emanating from q in the direction of v; it lies at a distance |v| = ((v, v)) 1 / 2 from q along this geodesic. 
(We use the notation ( , ) for the metric on A.) Thus, we find U q = exp Aq. We now modify our original definition of 
lA q . If the geodesies emanating from the constraint manifold C in the neighborhood of q flow to an arbitrary length, 
they will in general intersect one another. This can be witnessed even in the simple example of Sect. [□]. Thus, in 
defining 1A q , we assume that the geodesies flow for a small enough length to avoid such intersections and that this 
maximal length is independent of the point q on the constraint manifold. In summary, then, we foliate a neighborhood 
(which we call the tubular neighborhood) of the constraint manifold C by the transverse spaces U q , which we have 
in turn related to the normal spaces N q by the exponential map. Using the exponential map to construct tubular 
neighborhoods in this fashion is a standard technique. For details, see, for example, Lang |2(| and Vanhecke |2l| . 

Since we have identified normal vectors with points in the neighborhood of the constraint manifold, we view the 
transverse potential V± as a function defined on the normal spaces. With this interpretation, we will require that V±, 
as a function of q and the vectors in N q , be independent of q. By independent, we really mean independent modulo 
SO(d) rotations in Nq. Let us make this more precise. As in Sect. |nj, we specify the orientation of the transverse 
potential by an orthonormal basis E M , fj, = 1, d of the normal space N q . This basis forms a normal frame for the 
constraint manifold which we call the potential frame. For a given normal vector field u, we introduce the components 
u^, \x — l,...,d with respect to E M . The quantities coordinatize both the normal space N q and the transverse 
space U q , for which they are commonly called Riemannian normal coordinates [p2| . We use sans serif for the list of 
coordinates u = (u 1 , ...,u d ), reserving the bold notation u for the vector field. The neighborhood of C is therefore 
conveniently parameterized by (u,g). The heuristic constraint that V± be independent of position on the constraint 
manifold can now be made precise by the following statement: the transverse potential Vj_(u, q) as a function of (u, q) 
is required to be independent of q. 

In general, the construction of the parameters presented here is only possible locally on C. That is, it may be 
impossible to define in the neighborhood of the whole constraint manifold simultaneously. The construction can 
break down in two ways. First, it may be impossible to construct a tubular neighborhood for the entire constraint 
manifold. One can see this even with the simple example of Sect. [ll| If the one-dimensional curve spirals in on itself, 
then the width of the tubular neighborhood is forced to go to 0. (Recall that the width of the tubular neighborhood 
must be the same for all point on the constraint manifold.) Assuming that a tubular neighborhood does indeed exist 
for the manifold, there is still a second way in which the construction may break down. This occurs if there does 
not exist a potential frame E p which is globally defined. (This happens when the normal bundle is nontrivial.) For 
example, let the ambient space be a Mobius strip and let the constraint manifold be a curve which wraps around the 
Mobius strip once. Clearly, there does not exist a normal frame for C which is defined globally. It is our viewpoint 
that these two obstacles (in particular the first) are not common in physical problems. Even if one did encounter a 
problem in which the were not definable globally, since the results of this paper involve only the local form of the 
Hamiltonian, they would still apply to such problems. 



B. The Potential Twist Tensor 

In this section, we generalize the potential twist S, of Sect. ||, to a rank three potential twist tensor (also denoted 
S) defined for any q € C. For an arbitrary vector e E T q A, S e is a linear map on T q A. (Here, T q A is the (d + Tri- 
dimensional tangent space of A at q.) Let x £ T q A be an arbitrary vector tangent to C. Then, we define 

5 e x = 0. (4.1) 

Now let v 6 T q A be an arbitrary vector normal to C. We extend v to a vector field on C (defined in the neighborhood 
of q) by assuming that v is normal to C and furthermore that its components with respect to E M are constant. We 
now complete the definition of S e by prescribing 
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S e v = Pj_Vp lie v, 



(4.2) 



where V is the Levi-Civita connection p3|-|25[ on A and P± and P|| are the projection o pera tors onto the normal 
and tangent spaces of C respectively. It is straightforward to verify that S defined by Eqs. (hi) and (i.2) is indeed a 
tensor. 

Like the second fundamental form T (see Appendix 



(d, SJ ) 



S satisfies the antisymmetry property 
(f,S e d), 



(4.3) 



where d,e, f S T q A are arbitrary. To prove the above equation, we need only consider the case d = v G N q , 
f = w G N q , and e = x tangent to C, since all other cases vanish. Since S is a tensor, we may assume that v and 
w are vector fields and that their co mpon ents with respect to E M are constant. Since the frame E M is orthonormal, 
(v, w) is constant, and therefore Eq. (4.2) implies (v, S x w) = (v, V x w) = — (V x v,w) = — (w, 5 x v). 



C. Scaling of the Transverse Potential 

In Sect. ||, we imposed the constraint by taking the limit in which the width of the waveguide went to zero. Here, 
we describe a similar scaling procedure using, however, a more general transverse potential. Heuristically, we assume 
that V±(u;e) depends on the scaling parameter e in such a way that the potential grows narrower and deeper as e 
tends toward 0. To make this precise, we rescale the transverse potential in the following way 

V±(u;e) = e 2 V ± (eu;e), (4.4) 

where 

= en**. (4.5) 

We assume that has no dependence itself on e and that Vj_(u; e) is smooth in e at e — 0, by which we mean that 
Vj_(u; e) ca n be expanded as Vj_(u; e) = V^(u) +eVj_{u) + e 2 V 2 (u) + ... We also assume that V^(u) does not vanish. In 



Sect. VI B, we will make some very natural, further assumptions on the existence of bound states for the transverse 
potential and on the smoothness in e of the corresponding eigenenergies. 

As a concrete example take Vj_(u; e) = V±(u) to be a finite-depth square well with no e dependence. Then Vj_(u; e) 
is a finite-depth square well whose width scales as e and whose depth scales as 1/e 2 . Of course, these scaling rules 
apply to any potential satisfying the conditions described above. They guarantee that the expectation value of 
with respect to a transverse mode scales as e (assuming a fixed quantum number for the transverse mode). This fact 
shows that the wave function becomes more and more localized in the vicinity of the constraint manifold as e tends 
toward 0. 



V. EXPANSION OF THE KINETIC ENERGY 



The derivation of the constrained Hamiltonian (such as Eqs. (2.2) - (2.5)) proceeds in two steps. The first is to 
expand the kinetic energy in powers of e. The second is to transform to a basis of transverse modes and to apply a 
first order perturbation treatment to the expanded Hamiltonian. This section is devoted to the first step. 



A. Definition of the Vielbein 

We will express the kinetic energy in terms of a vielbein E , a = 1, d + m, on A. Appendix |c| gives the necessary 
background for this technique. To span the transverse dimensions, we take E M — d/du^, fi — l,...,d, where it is 
understood that, for the purpose of the partial derivative, the position q G C is held fixed. In selecting vector fields to 
span the remaining dimensions, we first choose an arbitrary set of orthonormal vector fields Ej, i = d + 1, d + TO, 
defined over C and tangent to C. We then use E M to Lie transport these vector fields into the tubular neighborhood 
of C. That is, we require the Lie derivatives with respect to E M to vanish, 

[E M ,Ei] = 0. (5.1) 



G 



We use the following notational scheme in this paper. The indices a, 6, c, ... range from 1, d + m and label the 
basis vectors E a and any components with respect to this basis. The indices fx, v, a, ... range from 1, d and label 
the vector fields E M = d/du^ and their related components. The indices i,j,k,... range from d+ 1, d + m and label 
the vector fields Ej and their related components. Except where otherwise noted, we employ the convention that an 
index a, 6, c, fj,, v, a, or k, ... is implicitly summed over when occurring twice in the same expression. 

For future reference, we present some facts regarding the structure constants 0^ b , defined by [E a ,E&] = 0^ b E lc - 
First, Eq. (5.1) immediately yields 



Furtherm ore, since E M = 
with Eq. (pj), we have 



Next we note that 



C ■ = C =0 

d/duV- is a coordinate basis on the transverse spaces U q , we find 0° 



bfi 



o. 



= [E M , [E j; ,E,]] = [E^-EJ = [d/Ou^^Ec 



(5.2) 

0. Combining this 
(5.3) 

(5.4) 



where the first equality follows from the Jacobi identity and Eq. (5.1) and the third equality follows from Eq. ([5.3]). 
We use the bracket notation [ ] in the final equality to emphasize that the differential operator acts only on the 
quantities inside the bracket. Since the E c form a basis, we have 



= 0. 



(5.5) 



Furthermore, since the Ej are tangent to C when restricted to C, we have 0fj = on C. Combining this fact with 
Eq. (5.5), we find that 



05 = o 



(5.6) 



within the tubular neighborhood. Thus, the distribution of vector fields Ej is integrable everywhere. (The submani- 
folds thus defined by Frobenius's Theorem are manifolds of constant V±.) 



B. Transformation of the Kinetic Energy 

As in Appendix the momentum operators are defined to be ir a — —iKE a . They are not in general Hermitian since 
the Hermitian conjugate is given by Eq. (C4). The kinetic energy is given by K — 7rjG afc 7r{ ) /2, where G a b = (E a ,E&) 
are the components of the metric tensor and G ab is the inverse of G a b- 

Appendix ^ also provides the framework for scaling the quantum wave function by an arbitrary (strictly) positive 
function s : A — > K (see Eq. (C5)) in order to modify the form of the kinetic energy. We apply this scaling formalism 
here, taking 



(5.7) 



where G — detG a fc. As mentioned in Appendix this scaling defines a new inner product of wave functions. We 
observe that the original inner product of two wave functions ip and ip' is given by 



Gv (p*(u,q)tp'(u,q), 



where v is the [d + m)-form 

v = E* 1 A ... A E< d+m *> = du 1 A ... A du d A E< d+l) A ... A E* {d+m \ 



(5.8) 



(5.9) 



Here, E* a is the basis of one-forms dual to the vielbein E a . We have also used the fact that E*^ = du^ . (Be careful 
not to confuse the notation ( | ) with ( , ) which denotes the Riemannian metric on A.) From Eq. (C6) we therefore 
find that the scaled inner product of two (scaled) wave functions ip and ip' is 



<v#'> 



v ip*(u,q)ip'(u,q). 



(5.10) 
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This scaled inner product gives rise to a scaled Hermitian conjugate, denoted f( s )- Using Eqs. (C4), (C8), and (| 
we find that tt^ is Hermitian with respect to the scaled Hermitian conjugate, 



Furthermore, Eqs. (C4), (C8), and (5.2) give the (scaled) Hermitian conjugate of Hj as 



tO) 



(5.11) 



(5.12) 



We now restrict the momentum operator itj to C. We write 7Tj|o to make this explicit; we use the notation |o for any 
quantity restricted to C since this corresponds to — 0. For present purposes, we consider the constraint manifold in 
its own right without being viewed as embedded in the ambient space. With this interpretation, the vector field ^j \o 
has a well-defined Hermitian conjugate which we denote by (^lo)^- This Hermitian conjugate is given by Eq. (|C4j), 
where it is understood that the symbol G now refers only to the determinant of the metric Gij on C. However, since 
the basis Ej|o is orthonormal, we have G = 1, and hence 



jk- 



(5.13) 



where the functions 7^ , defined on C , are the structure constants for Ej| . Since [E j 1 , Ej 1 0] = [Ej , Ej ■ ] 1 , the structure 
constants 7^ are equal to /3&|o< Comparing Eq. ( 5. IS ) to Eq. ( 5.12 ), we now have the following convenient description 



for 7r'- when restricted to C 



(-j (s) )L=(-iio) t - 



The scaled kinetic energy is given by Eq. (|C9|). Noting Eq. ( 5.11 ), we rewrite this as 

K »=\ (^G^TT, + TT^TTj + nl (B) G iV TT v + Tlj W G*^) + V s 

where, 



Vs = ~l (^G ab [n a lnG}[TT b \nG] + [ni^G ab [n b \nG}^j 



(5.14) 



(5.15) 



(5.16) 



We will henceforth drop the s index on K s , f (s), and ( | ) s , with the scaling being implicitly understood. 



C. Expansion of the Kinetic Energy 



In this section, we expand the kinetic energy Eq. (5.15) through order e°. Re call from Eq. (4.5) that is of order 



e , and hence the momentum 7r„ — ihd/du^ is of order e . From Eq. (5.1), we see that the mome ntum 7Tj is of 



order e°. Furthermore, from Eq. ( |5.5|) , we see that [3^ is of order e°, and combining this fact with Eq. ( 5.12 ) we find 

that 7rJ is of order e°. (Recall that the index "(s)" is now implicit.) These scaling properties imply that to expand 
Eq. ( 5.15 ) to order e°, we must expand V s , G y , G v , and G MI/ to orders e°, e°, e 1 , and e 2 respectively. 
Since E a is an orthonormal frame at — 0, we have the following identities 



G a b\o — G ab \a — 5 a b, 

Gab I s~i 
o- — — <**ab 



ib,<j\0, 



(5.17) 
(5.18) 



where we use the notation 
of G ab , 



a" for the derivative d/du a . Equations (5.17) and (5.18) yield the following expansions 



G ij {u) 
G«(u) 



= Si 



0(e), 



G^J u a + O(e 2 ) 
G^(u)=5 lw + G^ a \ u° + 



-G w> | w CT + O(e 2 ), 



-G-, 
1 



.\ u a u T + 0(e 3 ) 



5„ v - G^ a \ou a + -(-G^ + 2G tia , a G a ^ T )\ u rT u T + 0(e d ). 



(5.19) 
(5.20) 



(5.21) 
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The derivatives of the metric G a b appearing above are conveniently expressed in terms of the potential twist tensor 
S and the Riemannian curvature R on A. To see this, we first introduce the components of S, T, and R via 



Sabc = (E , Se^E;,), 
T a bc = (E a ,TE c Eb), 
Rabcd = (E a , RE c E d ~E>b), 

where T is included for completeness and for future reference. We then have the following identities 



G 



|o = 0, 



G 



1 



flV,(TT \0 



(5.22) 
(5.23) 
(5.24) 

(5.25) 
(5.26) 

(5.27) 



In actuality, the |o notation on S^crj is redun dant sinc e S i s only defined on C, but we will make use of this notati on a s 
a conve nient reminder. We will derive Eqs. ( 5.25| ) - ( 5.27 ) momentarily, but for now we insert them into Eqs. ( 5.19| ) 
- (5.21) to obtain 



G l >) 



0(e), 



G^(u) = -S MCTJ V ff + 0(e 2 ), 



G""(u) = ( 



S LLOk $V 



u a u T + 0(e 3 ), 



w here w e have use d the well-known symmetry of the Riemannian curvature R a bcd — Rcdab 
- (5.30) into Eq. ( 5.15| ) to arrive at the main result of this section, 



K = K \ 



K, 



yp 

v ex 



0(e), 



where 



= -7T 

2 



(5.28) 
(5.29) 

(5.30) 

We next insert Eqs. ( 5. 28] ) 
(5.31) 

(5.32) 
(5.33) 



K R = -R^ 
6 

V?=V s \n. 



loA^Ao 



(5.34) 
(5.35) 



We have taken advantage of the antisymmetry property Eq. (4.3), S^i = —S vll i, and the well-known antisymmetry 



relations R^vor — Ru^iar 



R^to to introduce the operators 



A„ 



1 



1 



(5.36) 



which are generalized angular momentum operators acting on the transverse space. That is, they generate SO(d) 
rotations i n the tra nsver se space. They are the generalization of the angular momentum A defined in Sect. [H| 

In Eqs. ( 5.32 ) - ( 5.34 ), we have employed the standard practice of raising tensor indices by contraction with G ab . 
Thus, 7T^ = G^na, S» vi = G^G^G 1 " S abc , etc. However, since G ab | = S ab , the raised components and lowered 
components of any tensor evaluated at = are actually equal. One could, therefore, equally well have written 
Eqs. (5.32) - ( |5.34 ) with all components lowered. The purpose of using raised components is simply to express these 
equations in manifestly covariant form. 

We now men tion a few facts concerning the Hermitian conjugate which we used to derive Eq. (5.33). First, notice 
from Eq. ( 5.11 ) that = A^. Also notice that since S^ vi \q has no dependence on u M , S^ v% \o and A^ v commute. 
This means in particular that the Hermitian conjugate in Eq. (5.32) may be applied to the 7r 4 |o term alone. Finally, 

we used Eq. ( 5.14 ) to relate (7Ti |o)^ to (7rj)|o. Notice from Eq. ( 5.15 ) that if Ei|o is a coordinate basis on C then the 

Hermitian conjugate may be dispensed with altogether; 

We call the terms K± , K?< and Kr appearing in Eq. ( |5.31 ) the transverse kinetic energy, the (preliminary) tangential 
kinetic energy, and the curvature energy respectively. The last term V^ x , being a scalar, nondiffcrential operator, we 
call the (preliminary) extrapotential. The three terms K^, Kr, and V£ x are all order e . The transverse kinetic energy 

K± is of order e~ 2 and therefore blows up as e shrinks to 0. The energy associated with this term will therefore be 
subtracted of f with the remaining three terms giving rise to the residual kinetic energy. Notice that each of the four 
terms in Eq. (5.31) is Hermitian with respect to the (scaled) Hermitian conjugate. 
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D. Proof of Identities (5^25) - (5.27) 



We return now to justify Eqs. ( 5.25| ) - (5.27). Considering Eq. (5.25), we have 



G 



(V Ect G m -)|o 
<V e „E m ,E j -)| 



(V Ect (E^E,))^ 
(E^V^E^L 



(5.37) 



where in the first equality, we replaced the coordinate derivative by the covariant derivative, treating G^j as a scalar 
function. The second equality is the definition of G^j, and the third equality follows from the Leibniz rule and the 
vanishing of the metric tensor under covariant differentiation. We next define the vector 



M CT/i = V E „E M , 



(5.38) 



which we now demonstrate vanishes on C. To prove this, we first show that it is everywhere symmetric in /i and a by 
using the general formula 



V d e- V e d= [d,e], 
where d and e are arbitrary vector fields on A. By substituting d = E^ and e 



(5.39) 



E M and recalling that [E^E,, 



0, 



we find that M„ 



M M(T . Since is symmetric, it vanishes if and only if /w^M^,, = for an arbitrary list of 

(constant) real numbers v a . For such an arbitrary list, we define the vector field v = v a E cr over A. Since dv^ /du a = 0, 
we see from Eq. ( |5.38 ) that v (J v^'M. ( j fi = V v v. Since the quantities u M are defined via geodesic flow away from C, 
an integral curve of v which passes through C is itself a geodesic. By the geodesic equation, (V v v)| = 0. Thus, 
-y CT - y ^M (TA1 |o = and hence 



(V Ect E m )| = 0. 



We return to Eq. (5.37) and write 



G w >|o = (E^, V E „Ej)| - (E M , V E jE ff )| =5 W -| , 



(5.40) 



(5.41) 



wher e th e first equality follows from Eq. (5.40), the second from Eqs. (5.3£) and (5.1), and the third from Eqs. (5.22) 
and Q. 

To prove Eqs. (15T261) and ( ^271) , we fix a point q on the constraint manifold and restrict our attention to a single 
transverse space U q which we temporarily forget is embedded in A. Recall that the vectors E^ are tangent to U q and 
G^u = (E M , 'Ely) is the metric tensor on U q . Furthermore, the coordinates are Riemannian normal coordinates on 
Uq and it is well-known that the expansion of the metric to second order in the Riemannian normal coordinates is p3] 



G^iu) 



V ~ -i? M <T l /r|o"' T W T 



(5.42) 



where R is the Riemannian curvature of the transverse space U q . The vanishing in Eq. (5.42) of the term linear in u 
proves Eq. ( |5. 26 ). Similarly, the quadratic term in Eq. (5.42) yields 



G^i/ : ar |o ^ (-^T^crz^T ~t~ Rvafir*) |o- 



(5.43) 



To complete the proof of Eq. ( |5.27 ), we must prove that the components R^ovt |o of the Riemannian curvature on 
U q agree with the components i? MCT i/r|o of the Riemannian curvature on A. To prove this, we use the Gauss relation 
given by Eq. (B10) and which we reexpress here in component form 



rp(L rp rpd rp 

1 crv-LaiiT ~ 1 aT 1 ajiv 



(5.44) 



Since we are applying the Gauss equation to the submanifold U q instead of C, we place an overbar on the symbols for 
the second fundamental form and the Riemannian curvature. Here, T is the second fundamental form oiU q . Recall 
that Py and P± were defined to be respectively the tangent and normal projection operators onto C. We extend the 
definition of these operators for m m not equal to by defining P|| and P± to be the normal and tangent projection 



operators respectively onto U q . With this definition, the second fundamental form T is given by (see Eq. (Bl)) 

TJ = PnV Pxe P±f + P x V P±e P\f, (5.45) 
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where e and f are arbitrary vector fields over U q which are tangent to A. Since E M is tangent to lA q everywhere, we 
have 



T ai iv = (E ct ,Te„E m ) — 0. 
Furthermore, since Ej is normal to U q at — 0, we have 

Ti pv \o — (Ej^E^Eju/lg = (Ej, Ve„E m )| = 0, 



(5.46) 



(5.47) 



where the last equality follows from Eq. (5.40). Combining Eqs. (5.46) and (5.47) yields T a J\q — T aflv \o — from 
which follows, using Eq. (5.44), -R MCT „ r |o = R^ V t\o- This concludes the proof of Eq. (5.27). 



E. Computation of the Extrapotential 



In this section, we analyze the extrapotential Vf x formed by evaluating Eq. (5.16) at = 0. As we will see, Vg x 
may be expressed solely in terms of the second fundamental form T of the constraint manifold and the Riemannian 
curvature R of A evaluated on C with no dependence on the potential twist S. Specifically, we will derive the following 
manifestly covariant form 



h 2 

V'P = - — 
8 



2TWT iw --T^T^. + 2i?% 



3 " 



(5.48) 



Setting Rabcd = 0, the above equation agree s wit h da Costa (Ref. jTj, Eq. (33)). Da Costa also assumes that 5 = 0. 
Since we do not make this assumption, Eq. ( |5.48|) is a generalization of da Costa's result to both the case of nonzero 
Riemannian curvature in the ambient space and nonzero twist of the potential. 

There are several other convenient forms for Vf x . We first introduce the following notation 

(5.49) 
(5.50) 
(5.51) 
(5.52) 

T^T^^T^T^, (5.54) 



■R 


nab 1 
- K ab\0, 


n± 








a 


= R ij ij , 


T 2 


^ rpabcrp 

— 2 a ^ C ~ 


M 2 


rpab rpC 

— 1 a 1 be — 



where we use Eqs. (pi[), (B6), and (B7) in Eqs. (5.53) and ( [3.54 ). The quantities 1Z and TZ are the scalar curvatures 
on A and C respectively. The quantity M is called the mean curvature. Using the fact that 1Z = TZ» +1Z± + 2i? v 



we rewrite Eq. (5.48) as 



if i. 



V p = 

ex 



2T 2 - M 2 



R-R\\- 



Furthermore, the Gauss Eq. (B1C) yields 



from which we find 



yp 



' 8 

rr- 



T = M -TZ + 1Z 



1 



11- -lis 



T 2 -K 

M 2 -2-R + TZ + TZi 



(5.55) 

(5.56) 

(5.57) 
(5.58) 



Assuming the tensor R = 0, Eq. ( |5.58 ) agrees with Ref. 0, Eq. (36) 



The remainder of this section is devoted to the derivation of Eq. (5.4S 
we observe that since G a b|o = &ab, we find 



Considering the first term of Eq. (5.16) 
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[ 7 r i lnG]| = 0, 

[Tr^lnGJIo = -ih(G ab G ba ^)\ Q = -ihG aa> „\ = -iWjj^lo, 
where we have used Eq. ( 5.26|) in the last step of Eq. ( 5.60 ). Equations ( |5.59 ) and ( 5.60 ) yield 

\ (G ab [7r Q lnG][^lnG]) | - ~ (G^G W , M ) | . 



(5.59) 
(5.60) 



(5.61) 



Considering the second term of Eq. ( |5.16| ), we note 

[7rtG afc [^lnG]]| = ( \x\ , G ab ] [n b In G] + G ab [ir^ [ir b In G]] ) | . (5.62) 

The first term of Eq. ( |5.62 ) vanishes from Eqs. ( 5.11 ), ( 5.12 ), ( 5.29 ), ( 5.3C ), and ( |5.59| ). The second term evaluates to 

(G ab [7rt[^lnG]]) | = [w^lnG}}\ = K 2 {G ab ^G ab4l - G 00lWt ) |„ 

2 



(5.63) 



where the first equality follows from Eqs. (5.11), ( 5.12 ), (5.17), and ( 5.59| ), the second equality is a straightforward 
computation, and the third equality results from Eqs. (5.2E]~^ ( [5.27 ). Collecting the preceding results, we find 



V p = -- 

v ex 



.Gu^Gjj^ 



G 



lt,fj,fj, 



(5.64) 



The various derivatives of Gij appearing in the above may be reexpressed using the following identities, to be derived 
momentarily, 



R 



Upon inserting Eqs. ( 5.65 ) and ( 5.66 ) into Eq. ( 5.64 ) one obtains Eq. ( 5.48j ) . 
Considering Eq. ( |5.65 ), it follows from 



G 



ij,fj.\0 



(Ve^E^E,) 



| = (Vb 4 E^Ej)|„- 



»MJ 1°' 



(5.65) 
(5.66) 



(5.67) 



wher e in the second equality we used the Leibniz rule and interchanged the derivatives by virtue of E qs. (5.39) an d 
( |5.l| ). The final equ ality follows from the definition of the second fundamental form Eq. (Bl) and Eqs. (B2) and ( |B5| ). 
Considering Eq. ( |5.66 ), we have 

Gij,nv\o = (VE„V E(J (E i! Ej)) | 

= (V E! E M ,V Ej E,)| Q + (Ve.E^V^E^)^ 

+ (V El/ V Ej E /J ,E,,-)| + (E i ,V El/ V Ei E /1 >| . (5.68) 

In the second equality, we again applied the Leibniz rule and interchanged derivatives by virtue of Eqs. ( [5.39] ) and 
(Of) . We next note that the covariant derivative of E M by Ej is given by 



(V Es E 



HJlQ 



(P||V Ei E M +P x V Ei E 



(T Ei E M + 



Wl0 ' 



(5.69) 



where the first equality follows from the fact that P\\ +P± is the identity and t he second from the definitions Eqs. (Bl) 
and ( p~2] ) and the fact that E M is normal to C. We also observe from Eq. ( 5.40j ) that (V Ei V El/ E M ) |o = 0, and therefore 

(V Ei ,V E ,E M )| = (i? EvEi E M )| , (5.70) 

wher e we have used Eqs. ( |5.l| ) and (|B^). Inserting Eqs. ( 5.69 ) and ( 5.7C ) into Eq. ( [3.68 ) yields the desired result 
Eq. (Of). 



VI. THE CONSTRAINED HAMILTONIAN 

In Sect. M we expanded the kinetic energy in e, obtaining two terms. One term, the transverse kinetic energy, is 
of order e , the other term is of order e°. In this section, we apply (degenerate) first order perturbation theory to 
derive a constrained Hamiltonian for the eigenenergies. In doing so, we introduce the transverse modes characterizing 
the wave function away from the constraint manifold. 
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A. Rescaling by e and the expansion of the Hamiltonian 



By adding the potential energy Vj_(u) to the kinetic energy Eq. (5.31), we have the following Hamiltonian 

H = H ± + H\\+0(e), (6.1) 

where 

H ± =K± + V±, (6.2) 

Hl = Kl + K R + V£ x , (6.3) 

are called the transverse and (preliminary) tangential Hamiltonians respectively. 

In order to clarify the subsequent perturbation analysis, we explicitly exhibit the e dependence of various quantities 
by rescaling them in e. To begin, we repeat the previous definition Eq. ([T^) of the rescaled quantities and also 
define rescaled momenta tt^, 

u>* = eu>*, (6.4) 
7T P = -fry,. (6.5) 

Notice that both and tt^ scale as e°. In general, the scaled version of a quantity (denoted with a tilde) is defined 
such that the lowest order nonvanishing term of its expansion in e is of order e . Thus, for a quantity homogeneous in 



e, the scaled version is independent of e. For convenience, we repeat the definition Eq. (4.4) of the rescaled potential 



energy V± and also define a rescaled transverse kinetic energy and transverse Hamiltonian 



V±(u;e) = e 2 V ± (u;e) = e 2 V ± (eu;e), (6.6) 
1 

2' 

H±(e) =e 2 H ± (e) = K± + V±(u;e). (6.8) 



K± = e 2 K ± = ^7f M 7r M , (6.7) 



By our previous assumptions in Sect. [V C| Vj_(u;e) is smooth in e and does not vanish at e = 0. As for K±, it is 



clearly independent of e. Thus, H±(e) is smooth in e at e = 0; its lowest order term is order e°, but depending on V±, 
it may have higher order terms as well. Recall that K$, Kr, Vf x , and ffj? are already independent of e and therefore 
need no further scaling. For notational continuity, however, we nevertheless define 

K\ = K\, (6.9) 

K R = K R , (6.10) 

Kl = V e p xl (6.11) 

H\ = HI (6.12) 

We rescale the full Hamiltonian by defining 

H(e) = e 2 H{e) = H x (e) + e 2 HP + 0(e 3 ). (6.13) 

In a typical Taylor series expansion of H(e), we would remove the order e term from H±(e) and leave it as a separate 
term. We would also combine the order e 2 term of H±(e) with the tangential Hamiltonian £ 2 H^. Here, however, we 

wish to keep the e and e 2 terms tog ether in the transverse Hamiltonian H±(e). We therefore define a new perturbation 



parameter k — e and rewrite Eq. (3.13) as 

H(e, K ) =H ± (e) + K HP + 0(e 3 ). (6.14) 



Our objective is to find the eigenvalues of H through order e 2 . Viewing e and k as formally independent in Eq. (3.14), 
our objective becomes finding the eigenvalues of H through second order in e and first order in k. Our procedure is to 
assume that the eigenvalues of H±(e) can be solved exactly (or at least through order e 2 ) and then apply first order 
perturbation theory in k. To simplify notation, we drop the e dependence (but not k dependence) for the duration of 
the derivation. 
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B. Transformation to the Transverse Modes 



The zeroth order term (in k) of H(n) is the transverse Hamiltonian H±, which has the form 

ti 2 d d 

Since H± depends only on the quantities u^, we may restrict its domain to functions of alone. For the moment 
we adopt this understanding for the domain of H± . We pick an eigenvalue E± (the transverse energy) of H± with 
finite multiplicity k and bounded eigenstates. We call these eigcnstat.es the transverse modes (with energy Ej_). We 
let Xn(u), n — 1, k, denote an orthonormal basis of these transverse modes. By orthonormal, we mean 

(Xn\Xn')u = I ' du 1 A ... A du d X *(u)xn'(u) = <W, (6.16) 



where the u subscript indicates integration only over the variables u 11 as opposed to the full (d+m)-form v in Eq. (|5.S|). 

We now adopt the understanding that H± acts on wave functions of both and q. Such an eigenfunction with 
eigenvalue E± has the general form 

k 

i;(u,q) = J2xn(u)Ml), (6-17) 

n=l 

where the functions 4> n (q) are arbitrary. We therefore identify an eigenfunction ip(u,q), having eigenvalue E±, with 
the k functions 4> n (q). Notice that with our current understanding for the domain of the operator H±, E± is a 
degenerate eigenvalue, even for k = 1. 

Recall the steps involved in first order degenerate perturbation theory. First, one either proves or assumes that 
the desired eigenvalue and eigenfunction is analytic in the perturbation parameter n. (Here, we simply assume this 
fact.) Next, one determines the zeroth order energy and zeroth order eigenstates. Then, one considers the operator 
formed by restricting the first order term of the Hamiltonian to the space of zeroth order eigenstates. The first order 
corrections to the energy are the eigenvalues of this restr icted operator. In the present case, the zeroth order energy 
is E±, and the zeroth order eigenstates are given by Eq. ( ]6. 17 ) . The first order correction to the Hamiltonian is kHZ. 



Denoting the first order correction to the energy by nEu , the eigenvalue equation for Eu is 



E 



Hu) cf> n ,= E«ct> n , (6.18) 

' nn' 



= 1 



where the (H^) nn r are the differential operators 

(#ll) nn , = (Xn \S%Xn' ) u = / du 1 A ... A du d X *n(u) (#fXn') (")• (6-19) 

We call (i?|| ) nn i the constrained, or tangential, Hamiltonian. 

We now recall that k = e 2 and reintroduce the explicit e dependence. Summarizing our analysis thus far, we have 
shown that an eigenvalue of H (e) through order e 2 is given by £ ± (e) + e 2 £||(e), where E±(e) and E\\{e) are eigenvalues 
of H±(e) and (H\\ )„„' (e). Of course, assuming smoothness in e, it is sufficient to solve for E±(e) and E\\(e) through 
orders e 2 and e° respectively We will therefore only require E^(e) and (-ff||) n „'(e) evaluated at e = 0, which we denote 
by and (H^) nn i respectively. Also, by virtue of Eq. (S.19), we assume for the remainder of the paper that the 



transverse mode s x n (u) are only order e° eigenfunctions of H±(e) 



We view Eq. (6.18) as a fc-dimensional vector wave equation for a vector wave function defined over the constraint 
manifold. We introduce the bold notation <fi(q) for the vector wave function with components 4> n (q) and the sans 
serif notation H|| for the matrix of differential operators with components (-ff||)rm'. Equation ( |6.18| ) can therefore be 
written more compactly as 

H [l (j) = E ll ^. (6.20) 

Havi ng c o mplete d the perturbation analysis, we have dropped the tildes from H|| and its eigenvalue Eu. Using 
Eqs. (6^), ( 5.33 ), ( 5.34 ) and a little algebra, we express H|| as 
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where 



H|| = K 1 1 + V ex , 



K|| = g^lol + S^'loA^J^TTil + S° T t \ a Kr), 



v ex = vn+ -s^s 



A (2) 

flVCTT 



QflVl Q(7T 



A(l T -A^A CTT ) + ii?^| Ag 



(6.21) 
(6.22) 

(6.23) 



and where I is the k x k identity matrix and and A^ I/0 . r are the k x k matrices having the following components 
respectively 

(A M „)„ n < = (x« lA^Xn' ) u , (6-24) 

(tflr)nn> = (Xn \^KrXn> > u ■ (6-25) 

Equations ( |6.21 ) - ( 6.23| ) encapsulate the main result of this paper. Specifically, the constrained Hamiltonian I — 1 1 1 is a 
k x k matrix of differential operators. It is the residual Hamiltonian remaining after the infinite transverse energy E± 
is subtracted off. The kinetic energy Km , which we call the (final) tangential kinetic energy, differs from the "standard" 
kinetic energy in that it has a gauge potential term. Physically, the gauge potential couples the tangential momenta 
to the generalized angular momentum of the transverse modes. The quantity V ' ex is a k x k matrix of nondiffcrential 
operators which we call the (final) extrapotential. Notice that any possible off-diagonal coupling in I — 1 1 1 is due to 
the angular momentum of the transverse modes. The preliminary tangential kinetic energy K?,, extrapotential Vf x , 

and tangential Hamiltonian are distinguished from the (final) tangential kinetic energy K||, extrapotential V ex , 
and tangential Hamiltonian I — 1 1 1 by the "p" superscript. We often drop the "preliminary" and "final" modifiers when 
referring to these terms, relying on their symbols and context to make our precise meaning clear. 

In the event of a nondegenerate transverse mode, that is k = 1, Hii becomes a scalar wave operator H\\ a cting on 



presented in Sect. VIII B 



scalar wave functions defined over the constraint manifold. In this case, we see the emergenc e of Eqs. (|2.3| ) - (|2.5| ) 
presented in Sect. 0. T he exact derivation of these equations from the more general Eqs. (6.21) - (6.23) will be 



C. Nonconstant Transverse Potentials 



Up to now, we have assumed that the transverse potential Vj_(u) is constant (modulo SO(d) rotations) along the 
constraint manifold C . For some physical systems this assumption holds exactly due to some symmetry on the ambient 
space, such as 5*0(3) rotations in the case of a rigid body. However, for other systems, this assumption may be only 
approximately satisfied; the constraining potential may in fact vary along the constraint manifold. This is true, for 
example, of a molecule evolving along a reaction path; there is no symmetry dictating that the frequencies of the small 
transverse vibrations be constant. The purpose of this section is to illustrate how small variations in the transverse 
potential may be easily included within our formalism. 

The key idea is to only allow dependence on q at order e 2 . Specifically, we assume the transverse potential can be 
expanded as 

V ± {u,q;e) = V^(u) + eVl(u) + e 2 Vl(u,q) + 0(e 3 ). (6.26) 



Applying this expansion to Eq. (6.8), an eigenvalue E± of H±_ (assuming analyticity in e) can be expanded as 

E x (q; e) = E° ± + eE\ + <?E\{q) + 0(e 3 ). (6.27) 

The first two terms of E± = E±/e 2 blow up as e goes to 0. However, these two terms are constant in q and may thus 
be subtracted off. The next order term Ej_(q) does depend on q and is of the same order in e as Hy. Thus, Ej_(q) 
may be combined with the extrapotential V ex in Hy to form the effective potential 

V e /(g) = y ex {q)+El(q)\. (6.28) 

This is the only modification which needs to be made to our formalism. Notice that the transverse modes x n (u) need 
not be modified since they are defined to be only order e° eigenfunctions of H± and hence are unaffected by the term 
El(q)- 
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VII. ANALYSIS OF CONNECTIONS 



Both the preliminary and the final tangential kinetic energies and K|| exhibit a gauge potential proportional to 
S a bc- In this section, we study the geometric origins of these gauge potentials and compute their curvatures. 

We begin by reviewing the connection on normal vector fields over C. We note that many equivalent definitions 
exist for the general concept of a connection. For the purposes of this paper, a connection is taken to be a covariant 
derivative operator which acts on some space of vector fields. For more background, see any of a number of standard 
references f23|-f25|. For the remainder of this section, v is an arbitrary normal vector field over C and x and y arc 
arbitrary tangent vector fields over C. The normal connection W N is defined by 

V x / v = P ± V x v = P ± V x P ± v, (7.1) 

where V is the Levi-Civita connection on A. Notice that v is itself a normal vector field. The curvature of V", 
denoted B N , is computed to be 

R N v — f\7 N \7 N \7 N \7 N \7 N 
#xy v - (, V x V y - V y V x - V [x>y] 

= P± (V x P ± V y - V y P ± V x - V [x , y] ) P ± v 

= P± (P xy - V X P,| V y + V y P|, V x ) P±v = P ± (P xy - T x T y + T y T x ) Pjv, (7.2) 



where the first equality is sim ply the definition of the curva ture , the second follows from Eq. (7.1), the third from 



noting Pi = / — P|| and Eq. (B8), and the forth from Eq. (Bl). As expected, the curvature depends only on the 
nature of the embedding of C (via the tensor T) and on the curvature of A. If we assume that the tensors B N and 
R vanish, then we obtain the class of embeddings considered by da Costa 0. For such embeddings, one can choose 
a po tent ial frame with vanishing twist, thus eliminating coupling between the transverse modes. (This follows from 



Eq. (7.5) below and the fact that for vanishing curvature, one can always find a frame for which the gauge potential 
vanishes.) Also, for a non- twisting potential frame, the submanifolds of constant potential are orthogonal to the 
transverse spaces lA q . Hence, at all points the restoring force is directed inward tangent to the lA q . 

It is instructive to compute the gauge potential explicitly for the connection V . For this computation we first 
choose an arbitrary orthonormal frame (not necessarily the potential frame) V^, fi = l,...,d, for each normal space 
N q . We denote the components of an arbitrary normal vector field v with respect to V M by v^. Then, the components 
of Vg.v are given by 

(V£v)" = E iV " + (A? r v v v , (7-3) 

where we have defined the gauge potential 

(Af)^ = (V^V&V*) = (V„,Ve 4 V„>. (7.4) 

Due to the orthonormality of the V M , {Af)^ is antisymmetric in ^ and v. The gauge potential can therefore be 
viewed as a one-form on C with values in th e Li e algebra so(d), which contains all antisymmetric d x d matrices. If 



we choose V p = E M , we recognize from Eq. ( [4.2[ ) that the gauge potential is related to the potential twist tensor by 

i^i )ia> — S^vi. (7-5) 

This result will be import below for analyzing Ktf and K||. 

We now consider a function ij)(y,q), such as the quantum wave function, defined in the neighborhood of C. (We 
use the bold notation v instead of sans serif used earlier because we wish to emphasize the dependence of ip on the 
normal vector and not on its components with respect to a given frame, such as the potential frame.) The connection 
V N which acts on normal vector fields gives rise to another connection V" p which acts on the function i/j(y,q). In 

order to define (V x p ?/>)(v, q), we first choose a path q'(a) such that q'(0) = q and (dqf /da)(0) = x. We then denote by 
v'(a) the unique normal vector at each point q'(a) satisfying v'(0) = v and (V^ dQ ,v')(a) = 0. Then, the connection 

V" p is defined by 



(vfV0(v,<z) = ^ 



</-(v' (a), </(«))■ (7.6) 



The transverse kinetic energy Kf, can be directly related to the covariant derivative V" p . To do this, it is useful 
to compute the gauge potential of V" p explicitly. As before, we consider an orthonormal frame and denote the 



1G 



components of v by v^. Then, the function ?/>(v, q) can also be interpreted as a function of (v, q), where v = (v 1 , v d ) 
is the collection of components. We therefore have 



(V«^)(v,,)=i 



a=0 



yj(v'(a),q'(a)) = — 

dtp 



a=0 



V>(v,</(a)) 



a=0 



= (Ey)(9)^r(v,9) + (E i V)(v ) g), 



(7.7) 



where in t he t hird equality, the derivatives d/dv^ and are understood to have q and v v held fixed respectively. 
From Eq. (|7j) and the condition V^.v' = 0, we find 



Inserting this result into Eq. (7.7) yields 

(Vj|£>)(v,g)= [(E l + 4 P )^\ (v,g), 

where 

Ap = (Af 

and where we have used the antisymmetry of (Af ) /J1/ to introduce the operator 



i 



'flu" 



_9_ 



(7. 



(7.9) 



(7.10) 



(7.11) 



Obviously if = E M , then A MI/ = —ihfl^. The relevance of V" p for Kf, is now clear. By choosing = E p and 



(7.12) 



applying Eqs. ([7J), ((fj), and ( |7.10| ) to Eq. ( |5.33| ), we see that 



Thus the preliminary tang ential ki neti c energy is just proportional to the Laplacian defined in terms of the connection 
VH p . (Compare Eq. ( \TT% to Eq. Q.) 

Considering Eq. ( [7.10[ ), we see that the two gauge potentials (Af r ) All/ and A| p differ only in their representation of 



so(d). For (Af ) 11U , we use a representation by (i x d antisymmetric matrices, whereas for Af p we use a representation 
by the operators f2 M „. Therefore, the curvature of the connections V w and V" p are also related by simply switching 
the representation of so(d). Hence the curvature -B" p of V" p is 



(7.13) 



We now consider a fc-dimensional vector- valued function 4>(q) with components 4> n {q)- The connection V" p induces 
a connection V" on (j> by the formula 



(Vj[ <f>)„ = ( X r, 



E 



(7.14) 



The tangential kinetic energy K|| is closely related to the connection V" as we now show. We take the orthonormal 
frame Vm t o be E^, a nd we recall that x„(u) is a function of alone and 4>(q) is a function of q alone. Then applying 
Eqs. Q, $7j|, and ( |7T0| ), we find 



v x An — Wxi "^ifArij 



(7.15) 
(7.16) 



where (5 X ) M „ = (E„,S X E„). Then Eq. (7.14) yields 
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Vjl = xl + All , 

A x = (Sx)'* ^/1J/J 



where ft^,, is the k x k matrix with components 

i 



(7.17) 
(7.18) 



(7.19) 



From Eq. (7.15), we note that the components of A x can also be written as 



(4) W = (Xn|Vib(n'>u. 



Equations (7.17) and (7.18) show that the tangential kinetic energy K||, Eq. (6.22), is given by 



(7.20) 



K|| = yKjMv 



(7.21) 



analogous to Eq. ( 7.12 ) for K?,. 

The connection V" is closely related to the adiabatic transport of quantum states and the associated geometric 
phase due to Berry p4fl . If a set of k degenerate quantum states £«(??), n = 1, depending smoothly on a set 
of to external parameters 77 = (771, f] m ), is subject to an adiabatic variation 7 7(a) of these parameters, then the 
^n(ct) = £n(f?(oO) satisfy (£ n \di; n i / da) = Yli(£n\d£n' /dr)i)(dr]i/da) = 0. Simon |26| recognized that this condition 
defines a connection 



V s 



_d_ 



l + Af 



(7.22) 



acting on the vector-valued wave function £ = (£ 1; ...,£&) parameterized by 77. The gauge potential A^ is a k x k 
matrix with components 



(Af)r 



drji 



(7.23) 



If the parameters r\i are themselves quantized, then the momentum conjugate to r\i is not simply —ih(d/dr)i)\, but 
rather —i^g/Q v . — —ih(d / dr)i\ + hf ). This situation applies, for example, to the Born-Oppenheimer theory of 
molecules, wherein the parameters m describe the positions of the nuclei and the £ n represent the quantum state 
of the electrons p7[ . For the constrained quantum systems considered in this paper, the ordering in e adiabatically 
separates the transverse modes \n (analogous to the £ n ) from the motion along the constraint manifold (analogous 



to the space of Therefore, the gauge potential A x occurring in Eq. ( [7.17| ) is essentially the same as Berry's 
gauge p otent ial Af occurring in Eq. (7.22). We say "essentially the same " bec ause the coordinate derivative d/drji 
of Eq. (7.23) has been replaced by the covariant derivative V" p of Eq. (7.20), this covariant derivative being the 
geometrically natural connection for the transverse modes. 

We next compute the curvature of the connection V". In terms of the gauge potential A" — (S' x ) A " y Q M ^, we have 

Bjly - (dA'l)(x,y) + [All, All] = (d5^)(x,y)Q^ + {S x )' M '(Sy) aT [Q tu/ , Q aT ], (7.24) 
where dS^ is the exterior derivative of , viewed as a one-form over C. We determine dS^ u from the formula 



Eq. (7.13) for the curvature B" p . We first note 

slip = (dAiif)(x, y ) + [4 p ,4 p ] = {dsnfoyWv + (s x r(s y r[fv^-L 



(7.25) 



where we have used Eqs. ([7^) and ( 7.10 ). It is straightforward to verify that the Q UL , satisfy the following commutation 
relations, 



(7.26) 



and hence Eq. (|7.25[) reduces to 



ft, is • 



(7.27) 



18 



Combining this equation with Eq. (7.13) produces 

(dS»")(xL, y) = (i? xy - T x T y + T y T x - S x S y + S y S x y. 



(7.28) 



Combining Eq. ( 7.28| ) in turn with Eq. ( 7.24 ), we arrive at the following useful formula for the curvature of V", 

bL = (flxy - T x T y + T y T x - s x s y + SyS^a^ + (s x y»(s y y T [q^, q ctt ]. (7.29) 



Using the commutation relations Eq. (7.26) the above equation can be recast as 



B\ y = (i? xy - T x T y + TyT^" <V 



(2) 

where Q-^o-t is the k x k matrix with components 



(^Lr)nn' = (Xn \^^MarXn') ■ 



(7.30) 



(7.31) 



VIII. SPECIFIC CASES AND EXAMPLES 



We consider several concrete examples to help clarify the general theory. 



A. Codimension One Case 



We assume here that the codimension of the constraint manifold is d — 1. Since there is only one normal dire ction , 
we expect the potential twist to vanish. Indeed, this follows from the antisymmetry property S^i = —S utl i (Eq. (fO|)) 
and the fact that p, = v = 1 . Similarly, the normal components of the Riemannian curvature R^o-t also vanish due 
to the well-known antisymmetry property R a bcd = —Rbacd = —Rabdc- From this fact follows TZ± = R H|o = and 
1Z = 72.ii. The expressions for T 2 and A4 2 can also be simplified by introducing the rank two symmetric tensor W 
defined on vectors tangent to C and with components W % ■ = T\ -. (This tensor is often called the Weingarten map.) 



Then T 2 = Tr (W 2 ) and M 2 = (Tr W) 2 . Hence, the tangential Hamiltonian Eq. (6.21) becomes 



h ii = K|| + y exi 

K|| = ^lo^l, 



y ex = v£j = --[T 2 -n + n n ) i 



[M - 211 + 272-u I 



(2T 2 - M 2 ) I, 



(8.1) 
(8.2) 

(8.3) 



where we have used Eqs. ( |5.55 ), (5.57), and ( [5.58 ). Notice that the tangential kinetic energy is proportional to the 
standard Laplacian on C. All reference to A Mt/ has vanished, and hence all coupling between the degenerate transverse 
modes has been eliminated. The /c-dimensional Schrodingcr equation therefore separates into k independent scalar 
Schrodinger equations. 

We consider the case where the ambient space A is a flat two-dimensional space and the constraint manifold C is a 
curve in that space. Then, we note that 72 = 72 = 72|| = 0. Furthermore, the second fundamental form, or equivalcntly 
the Weingarten map, has only one nonzero component. We denote this component by W = Wu — n = 1/p, where k 
is the external curvature and p is the radius of curvature. Then, the extrapotential is 



yp 



' 8 p 2 



(8.4) 



As in Sect. || the sign on Vf x is such that <fi is attracted to regions of high curvature. This extrapotential was derived 
earlier by Marcus ||] and Switkes, Russell, and Skinner 

We next consider the case where A is a flat three-dimensional space and C is a two-dimensional surface. We still 
have that 72 = 72|| = 0. Furthermore, the eigenvalues of the second rank two-dimensional tensor Wij are n\ = 1/ ' p\ 



19 



and K2 = l//?2) where p\ and p 2 are the two external radii of curvature. Then the extrapotential V^ x is conveniently 
written 



V* = -f [2Tr (H/ 2 ) - (Tr VF) 2 ] = - * (1 _ 1 
8 8 \pi P2 



■— (Ki - K 2 ) 



■5) 



This result was previously derived by Jensen and Koppe Q| as well as da Costa Q. 



B. Codimension Two Case 



We assume here that the codimension of the constraint manifold is d — 2. This allows us to define the quantities 
Si, A, and A< 2 ) by 



A,„, 



Ae„„, 



A (2) =A (2) 



where is the 2x2 antisymmetric tensor with £12 



-621 = 1- Furthermore, we have 



S.6) 
i.7) 



We express the tangential Hamiltonian Eq. ( J6.21 ) as 

h 1 1 = K|| + \/ ex , 

K,| = i( 7 r ! | l + 2S 4 | A)t(^l + 2S l A), 
Vex = V e p J + (2S t S l )\ (A {2) - A 2 ) + ^ ± |oA 2 . 



(8.9) 

(8.10) 
(8.11) 

(8.12) 



We consider the case of Sect. || where A is a flat three-dimensional space and C is a one-dimensional curve. First, 
we note 1Z = 7Z = TZ± = TZu = 0. Next, we denote the single component of tangential momentum by mi = %i. 

Since C is one dimensional, 7T|| = —ihd/da = nt, where a is the geodesic length. Furthermore, the potential twist 
is determined by the sole component S = Si. The second fundamental form can be identified with a normal vector 



T» = T** of magnitude k = l/p. Hence, M 2 = T^.T 3 1 , = = n 2 = 1/p 2 . Using Eq. (|5.58|), Eqs. (|Tlh and 



(3.12) therefore simplify to 



Km 



-(7r||l + 25A) 2 , 



V ex = ~^ K 2 l + (25 2 )(A( 2 )-A 2 ). 
8 

Assuming a single nondegenerate transverse mode, Eqs. and (|IJ) yield Eqs. Q and (|J). 



.13) 
.14) 



C. Rotationally Invariant Transverse Potential 



In this section, we assume the transverse potential Vj_(u) is rotationally invariant, depending only on the radius 
u = (it^it^) 1 / 2 in the normal space. The potential frame E M can therefore be any orthonormal frame we like. This 
freedom in the choice of potential frame produces a large range of possible potential twist tensors S, with the actual 
choice of S being simply a matter of convention. The Hamiltonian Hii in Eq. (6.21), however, should be independent 



(up to a rescaling of the wave function (p) of any such conventions. In the remainder of this section, we show explicitly 
how the dependence on S drops out of I — 1 1 1 under the assu mpt ion of rotational invariance. 
First, we observe that the transverse Hamiltonian Eq. (6.2) has the form 



H t = - 



1 d 



-i d 



A 2 



2 u^ 1 du 



.15) 
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where A 2 is the Casimir operator 



(8.16) 



Therefore, an eigenfunction Xn of H± is necessarily an eigenfunction of A 2 . We denote by Xn such an eigenfunction 

whose A 2 eigenvalue is A. A basic fact concerning the eigenspaces of the Casimir A 2 is that t hey b lock di agon alizc 

the generators A^. That is, (Xn\^~t*»Xn')u = if A 7^ A'.[ Based on the definitions Eqs. ( 6.24 ) and ( 6.25 ) for 

(2) . (2) 

A/jj/crr and A Mt/ , this fact implies that for the space of transverse modes for a given E±, njiva-T = A^A^,-. Similarly, 



^LiivaT = ft^fto-r. We therefore see from Eq. ( 6.23| ) that all S dependence drops out of V ex , 



V ex = VP x \ + -R^ aT \ /\^A aT . 



.17) 



Considering K , even though Eq. (6.22) is written in terms of the potential twist S, we showed in Sect. VII 
(specifically Eq. ( 7.21 )) that K|| can be expressed in terms of the Laplacian associated with the connection V". From 



Eq. (7.30) and the results above, we see that the curvature i?" of this connection is independent of 5 1 , 



B l y - ( R xy - T x T y + TyT x )' J,v Q^. 



(8.18) 



Now if two connections V" and V" have the same curvature, then their associated Laplacians can only differ by a 
rescaling of the wave function. Hence, the Hamiltonian H|| for different choices of the potential twist S can at most 
differ by such a rescaling. 



D. Harmonic Transverse Potentials 



We assume that the transverse potential is quadratic in the 

Vl(u; e )=]TiK( e ))VV< 



(8.19) 



and that the oscillation frequencies depend on e via oj^ie) = cD^/e 2 , with oj^ being independent of e. (For clarity, we 
make summation over the indices fj,, v, cr, ... explicit in this section.) We introduce the standard machinery of raising, 
lowering, and number operators for each degree of freedom, 

(8.20) 
(8.21) 




N,, = 



tl - 



Notice that and scale as e°. The transverse Hamiltonians H± and H± have the usual form 



(8.22) 

(8.23) 
(8.24) 



(8.25) 
(8.26) 



^his follows quickly from [A 2 , A M „] = 0. Note, (A - A')0c£|A M „X^)u = (Xn\ [A 2 , A M „]x^)u = 0. 
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and the transverse modes can therefore be labeled by the numb er o f quan ta n» in each degre e of freedom /i. We 
denote such a mode by Xn where n = (m, ...,rid)- Inserting Eqs. ( |B.2l| ) and ( [3.22 ) into Eq. ( 5.36 ) yields 



from which one quickly sees 

(Xn|A (ti „Xn) u = 
A significantly more involved computation yields 

h 2 

(Xn|A (tlI/ A CTT Xn)u = — g- 













-0 




-I) 



(8.27) 
(8.28) 
(8.29) 



We now assume that x = Xn is a n ondeg enera te transverse mode. The tangential Hamiltonian Eq. (6.21) is therefore 
a scalar operator. Using Eqs. ( |g.2£ ) and ( |3.29 ), H\\ is 

H\\ = K\\ + Vex, (8.30) 

(8.31) 



K \\ = 2 1T%][ ^' Kil 



h 2 ( l 

V eX — V^. — [ S^ % S llV i + ' 



1 -2 n 



.32) 



The most striking aspect of the above equations is that, due to the vanishing of (x|A^x)u, the tangential kinetic 
energy Ku is proportional to the standard Laplacian on the constraint manifold. Thus, all of the effects of external 
curvature and potential twist are contained in the extrapotential V ex . 



E. Potentials with Reflection Symmetry 



The vanishing of (x|A^i/X) u (and hence the potential twist as well) from K\\ in Eq. ( |S.31 ) follows from general 



considerations of reflection symmetry, and therefore occurs for a large class of symmetric potentials. 

Let Q S 0(d) be a reflection acting on the transverse coordinates u = (it 1 , ...,u d ), and assume that, for a given 
<7, u a is mapped to — u a and all other coordinates remain fixed. Thus, Q = Q _1 = Qt. Furthermore, assume that 
Vj_(u) is invariant under the action of Q, that is V±(Qu) = Kl(u). The reflection Q also has an induced action on the 
transverse modes, which we denote by Q and which is given by (Qx«)( u ) = Xn(Q -1 u)- Due to the symmetry of 
Q commutes with H±, 

[Q,H x ] = 0. (8.33) 

Furthermore, the following are easily verified 

Qif O} = (8.34) 
QttvQ 1 = -7r CT , (8.35) 
QA^Qt = -A CTM for all fi. (8.36) 



We now consider a single nondegenerate transverse mode, denoted simply by X- Due to Eq. (8.33), x m ust al so be 
an eigenfunction of Q with eigenvalue either +1 or —1 (since Q 2 = I). Combining these facts with Eq. ( 3.36 ) and 
recalling Q = Q\ we have 

<x|A CTM x) u = (Qx|A CTM Qx) u = (xn\QK^Xn) u = -(x|A CTM x) u , (8.37) 

and hence 

(x|A CTM x)u = for all [i. (8.38) 

If the potential Vj_(u) is symmetric with respect to at least d— 1 such Q reflections, possessing d — 1 distinct and 
orthogonal reflection axes then (x|A CT /iX)u vanishes for all /z, v = 1, d. For such highly symmetric potentials, K\\ 
is again given by Eq. ( |S.3l| ) and the only effect of the potential twist is to be found in V ex . This is the case for such 
common potentials as the simple harmonic oscillator, analyzed in the last section, as well as the d-dimcnsional square 
well. Note that this analysis says nothing about the off-diagonal terms of (A Mt/ ) rm / for a system with degenerate 
transverse modes; for such systems, there may indeed be a nonvanishing gauge potential. 
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IX. CONCLUSIONS 



We have rigorously derived the effective Hamiltonian of a constrained quantum system by considering the limit 
as the restoring force becomes infinite. In doing so, we have been careful to avoid unnecessary assumptions on the 
curvature of the ambient space, the form of the constraint manifold, and the manner of the constraini ng p otential. 



This general approach yields important new terms in the effective potential \/ ex , as outline d in Sects. VE and VI 
as well as a gauge potential in the tangential kinetic energy K||, as outlined in Sects. VI and VII. Furthermore, this 



general approach allows our theory to be applied to several examples of physical importance. These examples include 
reaction paths for molecular reaction and scattering problems, twisted quantum waveguides, the double pendulum, 
and models of polymers by rigid constraints. 

Perhaps the most important example of a constrained quantum system is the quantum rigid body. Though we lack 
space to include the analysis here, we have successfully applied our theory to this case. Physically, we have in mind such 
systems as semirigid molecules. If we assume that the standard Born-Oppenheimer ordering for semirigid molecules 
is valid, then our constrained Hamiltonian reproduces (through the lowest three orders in the Born-Oppenheimer 
ordering parameter) the standard results for the rotation-vibration energy levels of a semirigid molecule. (See, for 
example, Papousek and Aliev pq|.) For such molecules the gauge potential term in K\\ vani shes du e to the harmonic 
form of the constraining potential. (We assume a nondegenerate vibrational state; see Sect. |VniPj| .) For this reason, 
a more interesting example would be one in which the standard semirigid analysis breaks down. This occurs, for 
example, in rigid clusters of molecules held together by van der Waals forces. For these systems, the gauge potential 
will not in general disappear and should have measurable effects on the rotation-vibration spectrum. We will pursue 
these issues in future publications. 
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APPENDIX A: A BRIEF REVIEW OF CURVES IN R 3 



We cite a few important facts about curves in R 3 which we need in the body of the paper. For greater depth, see, for 
example, Spivak p3 . Consider a curve x(a) in R 3 . The parameterization of the curve is given by a which measures 
the arclength along the curve. Hence the tangent vector t = dx/da is of unit length. We denote the principal normal 
and the binormal by n and b, respectively. They are given by 



dt/da 
\di/da\ 
txn. 



(Al) 
(A2) 



The vectors (t, n, b) form an orthonormal right-handed frame. The derivatives of this frame are given by the famous 
Serret-Frenet formulas which may be summarized as 



d 

da 



' t ' 




K 0" 




" t " 


n 




-K T 




n 


b 




-T 




b 



(A3) 



where K,{a) and t{q) are called the curvature and torsion respectively. The curvature and torsion have units of 
reciprocal length. The reciprocal of n is the radius of curvature p — 
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APPENDIX B: THE SECOND FUNDAMENTAL FORM 



The external curvature of a submanifold C embedded in a manifold A is conveniently specified by a rank three 
tensor T called the second fundamental form. Since the second fundamental form is of critical importance in the body 
of this paper, we briefly review a few of its relevant properties. For greater detail, see Refs. ^9| , p0| . 

Throughout this appendix, d, e, f denote arbitrary vector fields tangent to A and defined over C; w, x, y, z denote 
vector fields tangent to C; and v denotes a vector field normal to C. The second fundamental form applied to e and 
f , denoted T e f , is a vector field defined by 

TJ = P ± V P||e P||f + P||Vp ||e P±f, (Bl) 

where V is the Levi-Civita connection on A and P|| and Pj_ are respectively the tangent and normal projection 
operators of C It is straightforward to verify that T is in fact a tensor. Furthermore, the second fundamental form 
satisfies the identities 

(d,T e f)=-<f,T e d>, (B2) 
T x y = T y x, (B3) 

wher e ( , ) denotes the Riemannian metric on A. In terms of the components T a i, c = (E a , Te c E;,) introduced in 
Sect. VC, we have 



^abc 


Tbaci 


T 


— T 
1 aji i 


^abfi 


= o, 


T 


— T — 



(B4) 
(B5) 
(B6) 
(B7) 



wher e th e first two equations are simply component forms for Eqs. (B2) and (JB3|) and the last two follow easily from 



Eq. (Bl) 



In Sects. VD and VE, we need the Gauss equation, a well-known identity relating the second fundamental form T, 



the Riemannian curvature R of C, and the Riemannian curvature R of A. The Riemannian curvatures are defined by 

Pdef - (V d V e - V e V d - V [d . e] )f , (B8) 

P xy z= (V x V y - V y V x - V [x , y] )z, (B9) 
where V denotes the Levi-Civita connection on C. The Gauss equation is then [p9| 

(w, P xy z) = (w, P xy z) + (T x z, T y w) - (T y z,T x w). (BIO) 

APPENDIX C: THE QUANTUM KINETIC ENERGY WITH RESPECT TO A VIELBEIN 

We present two expressions for the kinetic energy of a quantum system on a Riemannian manifold of dimension n. 
These expressions differ in the scaling of the quantum wave function. We refer the reader to earlier related analyses 
pi] -|33| for derivations and discussion. 

We express the kinetic energy in terms of a vielbein. By a vielbein on a Riemannian manifold, we mean a set of 
vector fields E a , a = 1, ...,n, forming a basis of each tangent space. The structure constants fi c ah of the vielbein are 
defined by 

[E 0) E 6 ]=^ 6 E C , (CI) 

where [ , ] denotes the Lie bracket. The structure constants vanish if and only if the vielbein is a coordinate basis, 
that is if and only if there exists a set of coordinates x a such that E a = d/dx a . We denote the components of the 
Riemannian metric with respect to the vielbein by G a b and the inverse of G a b by G ab . 



2 Our definition of the second fundamental form differs in the choice of domain and range from that in Ref. [Q. We follow 
the definition of Ref. H. 
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We define the kinetic energy of the quantum system by K = —h 2 A/2, where A is the Laplacian. In terms of the 
vielbein, the kinetic energy is J3l| , 



2 



where 



-i?iE n 



(C2) 



(C3) 



are the momentum operators. In the above f denotes the Hcrmitian conjugate. In general, the momenta ir a are not 
Hermitian. They do, however, satisfy the following useful identity 



7rJ = 7T a 



7r a In VG 



(C4) 



where G = detG a h. The bracket notation in Eq. (C4) indicates that the quantity inside the brackets is a scalar; that 
is, n a acts only on In \/G. 

Often it is useful to scale the original wave function ip by some real positive function s to form a new wave function 

<l>, 



tp = sip. 



(C5) 



Such a scaling produces a new kinetic energy operator acting on the new wave function tp. By conveniently choosing 
the scale factor s, the new kinetic energy may acquire a more convenient form than the old kinetic energy. To 
demonstrate how the kinetic energy transforms, we first observe that the scaled wave functions have a different inner 
product than the unsealed wave functions. Denoting the unsealed inner product by ( | ), the scaled inner product 
( I ) s is defined by 



1 



(C6) 



for arbitrary wave functions ip an d "0'- This scaled inner product in turn defines a scaled Hermitian conjugate 
of an operator A. Specifically, 



(C7) 



Applying Eq. (C7) to Eq. (C4), we find 



r t0) 



7I-T - 2[7r a lns]. 



(C8) 



The scaling of the wave function transforms the kinetic energy operator K into K s = sK(l/s). It can be shown 
5l| that K* reduces to 



K, 



1 



Trt^G ab n b + V s , 



(C9) 



where 



V, = ~ (G ab [n a lns][7r b lns] + [^0^ Ins}] 
= I (G Qb [7r a ln S ][^ln S ] - [7rtG a6 [7r 6 lns]]) . 



(CIO) 



FIG. 1. A twisted quantum waveguide. The cross sectional shape of the tube is constant and is chosen to be a triangle with 
no reflection symmetry. (Reflection symmetry would force (A) to vanish.) The vectors Ei and E2 determine the orientation of 
the sides of the waveguide, and a measures the distance along the axis. 
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